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Abstract 

We show that the existence of a pair of zero-energy modes bound to a vortex 
carrying a 7r-flux is a generic feature of the topologically non-trivial phase 
of the M — B model, which was introduced to describe the topological band 
insulator in HgTe quantum wells. We explicitly find the form of the zero- 
energy states of the corresponding Dirac equation, which contains a novel 
momentum-dependent mass term and describes a generic topological transi- 
tion in a band insulator. The obtained modes are exponentially localized in 
the vortex-core, with the dependence of characteristic length on the param- 
eters of the model matching the dependence extracted from a lattice version 
of the model. We consider in full generality the short- distance regularization 
of the vector potential of the vortex, and show that a particular choice yields 
the modes localized and simultaneously regular at the origin. Finally, we 
also discuss a realization of two-dimensional spin-charge separation through 

the vortex zero-modes. 
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1. Introduction 

Topological band insulators (TBIs) have recently opened a new frontier 
in both theoretical and experimental condensed-matter physics due to their 
peculiar properties (extensive reviews are in Refs. [HE])- These stem from 
the fundamental fact that TBIs are described by topological field theories 
[3], thereby taking the interest in studying them far outside the standard 
condensed-matter domain. The novelty of TBIs is their protection by time- 
reversal symmetry (TRS), which leads to a non-trivial topological Z2 invari- 
ant |4j of their free electron crystalline band structure, and to a description 
in terms of the topological BF theory in (2+l)D and (3+l)D [5]. Striking 
consequences of their topological nature are the remarkable effective field 
theories describing the responses of TBIs, for instance, axion electrodynam- 
ics in (3+l)D [0] and gravitational Chern-Simons in thermal response [HIE]- 
Such theories are governed by anomalies, and the possibility of their direct 
study, e.g. through the Witten effect of axion electrodynamics [8], is of great 
interest and potential for both the high-energy and condensed matter com- 
munities. Further beyond the non-interacting case, there has been a proposal 
of (3+l)D fractional TBIs [9J, having TRS and an axion angle different from 
or 7r; these have a description in the form of deconfined non-Abelian gauge 
fields explicitly realized using holography [10J. 

The salient feature of TBIs, crucial for their characterization and detec- 
tion, is that they are fully gapped in the bulk while possessing on their bound- 
ary gapless propagating modes protected by TRS. [H|2] The presence of TRS 
however limits, fundamentally and especially experimentally, the availabil- 
ity of robust probes of such bulk-boundary correspondence, aforementioned 
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anomalous responses, and Z 2 topological order itself. For instance, in (2+l)D 
TBIs, on which we are focusing from now on, the charge Hall response van- 
ishes, and instead a much more involved TRS invariant quantum spin Hall 
(QSH) effect characterizes the topological phase. It has been understood 
through numeric studies that a 7r-fiux vortex, which actually preserves TRS, 
can play exactly the role of a Z 2 probe in a QSH insulator through appearance 
of topologically protected zero-modes. [51 [EH H21 US] 

In this paper, we analytically study the properties of 7r-ffux vortex in 
presence and absence of Z 2 order. The general way by which we achieve this 
is using the M — B model initially constructed to describe the HgTe quan- 
tum well QSH insulator. [J?l [T5] The salient and universal feature of the 
low-energy (continuum) version of this model is that it describes a topolog- 
ical phase transition between a trivial and non-trivial Z 2 topological phase, 
through a massive Dirac-Schrodinger theory. This field theory, especially in 
the presence of a U(l) vortex, has not been widely studied for its own sake. 
A peculiar property of this theory is that the presence of both linear and 
quadratic kinetic terms together with the ordinary Dirac mass term allows 
for a gap-closing transition which changes the Chern number of the bands 
and the Z 2 invariant. The same theory turns out to harbor analytically 
solvable zero-modes tied to 7r-flux vortex, but only in the non-trivial phase. 

The relationship of the 7r-flux modes to the QSH phase and the ques- 
tion of their protection are general problems in the context of zero-energy 
fermionic modes bound to a topological defect. Namely, as Aharonov and 
Casher have shown in Ref. [16J, when non-relativistic (Schrodinger) fermions 
are coupled to a magnetic flux carrying n flux quanta there are precisely 
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n zero-energy modes in the spectrum of the Hamiltonian. Later, Jackiw in 
Ref. p2] has demonstrated that a magnetic flux with n flux quanta hosts 
precisely n zero-energy modes when coupled to relativistic Dirac fermions, 
and they are related to an index theorem for the Dirac Hamiltonian defined 
on a compact space. [H] On the other hand, the existence of fermionic zero- 
modes bound to a vortex in the complex scalar order parameter has actually 
been established in one and two spatial dimensions in the pioneering works 
by Jackiw and RebbifTU] and Jackiw and Rossi [20]. respectively. Their ex- 
istence is, at the deep mathematical level, tied to an index theorem that 
relates the spectral asymmetry of the corresponding Hamiltonian defined on 
an open space and a topological invariant of the background scalar fields. [21] 
At the same time, the results of Jackiw and Rebbi have been applied to the 
polyacetilene system through the Su-Schrieffer-Heeger model. [22] However, 
the connection between Jackiw-Rossi and the M — B model, relevant for the 
quantum spin Hall system, has been shown only very recently in Ref. [23J. 
Namely, these authors have demonstrated that the Jackiw-Rossi model in 
the limit when Zeeman coupling and the chemical potential are large reduces 
to the M — B model, and therefore the results presented in this work are also 
relevant for this problem. Moreover, a Hamiltonian of the same form as the 
M — B Hamiltonian describes non-relativistic p x + ip y superconductor, and 
the results of our work are thus relevant for this system as well. In partic- 
ular, our solution for the zero-energy mode bound to a 7r-flux vortex in the 
quantum spin Hall state corresponds to the Majorana state in the core of a 
7r-flux vortex in this topological superconductor. 

The flux-carrying vortex is by its definition singular in real space. It is well 
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known that a Hamiltonian with singular potential (e.g. Aharonov-Bohm flux 
vortex [MIES], Coulomb potential [261E3EE], delta function potential [29]), 
once made Hermitian through a self-adjoint extension [30], can exhibit finite 
or zero-energy bound states, even if the original Hamiltonian was scale-free. 
Therefore, the question of regularizing the vortex singularity, and thereby 
completely defining a Hermitian fermionic theory, becomes physically rel- 
evant. It fixes the real space profile of zero-modes at the vortex and the 
scattering phase shift there. This problem has been considered for (2+l)D 
[2HET], and (3+l)D [32] Dirac particles, as well as quasiparticles in super- 
conductors [25J . The Dirac-Schrodinger theory of the M — B model has not 
been studied in this context before, and we find here the general form of the 
regularized theory in presence of a magnetic 7r flux vortex. Furthermore, we 
show that a regularization corresponding to a thin solenoid limit surprisingly 
leads to localized zero-modes that vanish at the origin. 

In this paper, we first explicitly introduce the tight-binding and contin- 
uum M — B models. Then we use the continuum M — B model to ana- 
lytically show that the 7r-flux vortex hosts precisely a pair of exponentially 
localized zero-energy modes, and therefore the states found numerically in 
Refs. [33l [13] are indeed a generic feature of the M — B model. Moreover, 
the relationship of these midgap states with the topological properties of the 
quantum spin Hall state are also considered. Namely, we show that these 
modes, in fact, exist in the entire range of parameters describing topologi- 
cally non-trivial phase in the M — B model with the gap opening near the 
zero momentum (T point) in the Brillouin zone. For obtaining the explicit 
form of the midgap states, a short-distance regularization of the Hamiltonian 
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is necessary due to the singularity of the vortex potential. We show that a 
particular choice, corresponding to a limit of infinitely thin solenoid, yields 
the modes non-singular at the origin. We then consider in full generality the 
problem of consistent regularization of the vortex potential, and show that 
the corresponding self-adjoint extensions are given in terms of parameters 
described by U(2) matrices. Next, we test the prediction for the localization 
length of the zero-modes obtained from the continuum theory as a function 
of the transition driving parameter by comparing it to the results from a 
lattice regularization of the model. Finally, we discuss the quantum numbers 
of the obtained zero-energy states. 

2. The M — B tight-binding model on the square lattice 

We begin by considering a minimal tight-binding model proposed to de- 
scribe a two-band quantum spin Hall insulator [13] 

H = 5>t(k) I ^ (k) U(k) (i) 
V V o #*(-k) ) 

where \1/ T = (u^, v^, u±, v±) = (^-f, with u and v representing two low- 
energy orbitals. The upper and the lower blocks in the Hamiltonian are 
related by time-reversal symmetry, and H(k.) acting in the orbital space has 
the form 

H{k) = <vk(k), (2) 

where <r M , fi = 1,2,3, are the Pauli matrices, d li2 (k) = Asmk xy , and d 3 = 
M— 25(2— cos k x — cos k y ), the length unit is set by lattice constant a = 1, and 
summation over repeated indices is assumed hereafter. We also set h = c = 
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e = 1 in the following, unless otherwise stated. Since the above Hamiltonian 
has spectrum E(k) = ^Jd^d^ doubly degenerate in spin space, the band gap 
closes at the T-point (k = 0) in the Brillouin zone when the value of model 
parameters is M/B = 0. We will only consider the range of parameters 
< Mj B < 4 in what follows, so that the spectrum is gapped. This lattice 
model then describes a topologically nontrivial state with a Kramers' pair of 
counterpropagating modes on the edge of the system leading to a quantized 
spin Hall conductance [T5] a^ y = 2^-. For negative values of M/B, the model 
describes a trivial insulator (cr~ = 0). 

3. The Dirac-Schrodinger continuum theory 

By taking the continuum, i.e. large wavelength limit (|k| <C 1) of ([2]), we 
arrive at a Dirac Hamiltonian which besides the ordinary Dirac mass term 
(M) contains a Schrodinger kinetic term (B) 

H cS (k) = i 7o7i fcj + (M — £?k 2 ) 7o , (3) 

where the four-dimensional 7 - matrices are given by 7 o = 03 <8> t , 7l = (x 2 <8> 
r 3 , and 72 = —o\®Tq. Here, Pauli matrices {r , r M } act in spin space, 
with (To, to as the 2x2 identity matrices. The 7 -matrices satisfy canonical 
anticommutation relations { 7cn 7( g} = 25 a/ 3, with a, /3 — 0, 1, 2. Notice that it 
is enough to focus on a single spin projection, since the two spin projections 
are related by the time-reversal operator T = —1T2K, with K as the complex 
conjugation. For convenience, we have divided through Eq. ([T]) by the energy 
scale A (the lattice intersite hopping energy) and by the length-scale a, before 
redefining M/(Aa) — > M, Ba/A — > B, so that the continuum theory Eq. (J3]) 
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has a dimensionless parameter MB and parameters B and yVB/M with 
dimension of length (we revert to lattice units for comparison to the tight 
binding model in Section [7]). 



4. Zero-energy states 

Let us now consider the effect of the magnetic 7r-flux inserted into the 
system. As usually, the vector potential is coupled to the electronic degrees of 
freedom through the minimal substitution, k — > k + A, and the Hamiltonian 
^ for spin up electrons assumes the form 

H eS (k, A) = a^h + A { ) + [M — B(k + A) 2 ]a 3 . (4) 

The vector potential 

A= -,e. + xe y (5) 

represents the magnetic vortex carrying the flux $ = tt. Notice that the spin 
down electrons are coupled to the 7r-flux with the opposite sign because of 
the time-reversal symmetry, and in that respect 7r-flux thus acts on the spin 
components as a pseudomagnetic vortex in graphene does on the two valley 
degrees of freedom. [3H EH EE] Of course, the time-reversal invariance of the 
Hamiltonian Q is present only when the flux corresponding to the vector 
potential A is equal to 7r or 0. 

We now show that the above Hamiltonian possesses precisely one bulk 
zero-energy state with spin up. Time-reversal symmetry then implies the 
existence of the zero-energy state for electrons with spin down. Expressing 
the Hamiltonian Q in polar coordinates (r, cp), taking into account that 
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A = (l/2r)e v , we obtain 
H cS = -ie 

+ 



-up 



d r - -Bu 



a + — ie 

1 X 



lip 



M + B [ di + -d r + -di 



d r + -cL 
r 



0"3, 



<7- 



(6) 



where d v = dp + i/2, and a± = (o"i ± ia 2 )/2. It is easy to see that in case of 
an arbitrary flux $, the Hamiltonian Q also acquires the form (|6]), but with 
the operator d v = d v + i$/27r. 

In the presence of a vortex carrying a 7r-flux, we seek the zero-energy 
modes of the form 



e iltp vi(r) 



(7) 



where / G Z is the angular momentum quantum number, and the functions 
u, v are the solutions of the following equations 



A,_iii(-i(r) -i ( d r + 



I i Or - " ^ 



■uj(r) = 



tti-i(r) + A z+ i^(r) = 0. 



Here the operator A; is defined as 



A, = M + B [ d 2 r + -d r - ^ ) = M + BOi. 



Acting on Eq. (8) with the operator A i+ i, and using the identity 

[A I A + ;] = -(2/-l)^(ft. + ;) ) 

we can eliminate the function vi(r) from the same equation to obtain 

2BI 

A I+ iA,_i - 0,1 + -j-A,_i ) u^r) = 0. 



(8) 
(9) 

(10) 



(11) 



(12) 



After some algebra, the above equation may be rewritten as 



M 2 + {2MB -1)0^1 + B 2 Of_i u,_ 1 (r) = 0. (13) 

This result may also be obtained by noting that if the spinor in Eq. ([7]) is an 
eigenstate with the zero eigenvalue of the Hamiltonian Q, then it is also an 
eigenstate with the same eigenvalue of the square of this Hamiltonian. Using 
Eq. Q, one then readily obtains 

# e ff(k, A) 2 = B 2 (k 2 ) 2 + (1 - 2MB)k 2 + M 2 , (14) 

with k = k + A, and the operator k 2 after acting on the angular part of 
the upper component of the spinor (J7J) yields Eq. (13). Similarly, it may 



be shown that the function vi(r) in the spinor given by Eq. ([7]) obeys an 



equation of the same form as (13) with I — >■ l + l. From Eq. ( JX3[ ) we conclude 
that the function u/_i(r) is an eigenfunction of the operator 0;_i/2 with a 
positive eigenvalue 

0, 4llH (r)^Vi(r), (15) 

since the operator k 2 when acting on a function with the angular momentum 
I is equal to — 0;+i/2, and the eigenstates of the operator k 2 with a negative 



eigenvalue are localized. Eqs. (13) and (15) then imply 



2B K J 

and the function ui(r) ~ I t _ i (Ar) with Ii(x) as the modified Bessel function of 
the first kind. However, from the above solutions the only square-integrable 
ones are in the zero angular-momentum channel since Ii(x) ~ as x — > 0. 
Furthermore, for / = only the linear combination I\/2{x) — I-i/2{%) ~ 
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x~ 1 l 2 e~ x has the asymptotic behavior at infinity consistent with a finite norm 
of the state. In the above equation we should distinguish two regimes of 
parameters, < MB < 1/4 and MB > 1/4, for which the argument of the 
square-root is positive and negative, respectively. 

For < MB < 1/4, since the argument of the square- root in the above 
equation is positive, we obtain two zero-energy solutions 



*±(r) 



,-X±r 



2ttA±V 



(17) 



and, of course, A± > because of the square-integrability. On the other 
hand, when MB > 1/4, up to a normalization constant, the solutions have 
the form 



*i(r) 



* 2 (r) 



e -V?«»*cos (ryj§ cos 6 s ) 



e -rVfcos0 sin (r^McosQ^ 



-up 



(18) 



where 



However, since the identity 



1 J\\ -4MB 
- arctan - - 

2 1 - 2MB 



1 J\l-4x\\ 1 
x cos I - arctan — = - 

2 l-2x I 2 



(19) 



(20) 



holds for 1/4 < x < 4, the localization length of the zero-modes for MB > 
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1/4 is actually independent of M, namely Eqs. (18) become 

e-5Bcos(^) / e 



-up 



*i(r) 



i 



e~^B sin (^) / e-^ . 
tf 2 (r) = . (21) 



r 



/ 



Therefore, we can conclude that the Hamiltonian Q possesses zero-energy 
modes in the entire range of parameters M and B for which the system 
is in the topologically non-trivially phase, < M/B < 4. In particular, 



as it can be seen from Eq. (16), when 4MB < 1 zero-energy states are 
purely exponentially localized, while for 4MB > 1 the exponentially localized 
solutions have an oscillatory part with a characteristic length-scale exactly 
equal to the localization length. 

Notice also that in the regime when < MB < 1/4, there are two charac- 
teristic length scales associated with the midgap modes, £± ~ A^ 1 . Of course, 
after a short-distance regularization is imposed, only a linear combination of 
the two states survives. The physical interpretation of the two length scales 
depends on the form of the superposition of the state after the regularization 



has been imposed, as it may be easily seen from the form of the states (17). 
In the regime MB > 1/4, the zero-energy states are characterized by a single 
length-scale £i oc ~ 2B, which is at the same time the localization length and 
characterizes the oscillations of the exponentially decaying state. 

Therefore, the appearance of the zero-energy states bound to a 7r-flux 
vortex is a generic feature of the Hamiltonian ^ describing the quantum 
spin Hall system. Furthermore, in the vortex-free system, it may be shown 
by imposing open boundary conditions on the wave-function at one of the 
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edges of the system, for instance the one perpendicular to the x-axis, and 
at infinity, ^(x,y = 0) = *&(x,y = oo) = 0, that the Hamiltonian gives rise 
to gapless edge modes with the penetration depth given by exactly the same 
expression as the localization length for the zero-energy modes bound to the 
7r-flux vortex. The bulk-boundary correspondence may be thus probed by 
inserting a 7r-flux vortex in the quantum spin Hall system. 

5. Thin solenoid regularization of vortex 



The zero-energy modes, given by Eqs. (17) and (21 ), form an overcomplete 
basis in the zero angular-momentum channel, because the Hamiltonian Q is 
not self-adjoint, which is due to the singularity of the vortex vector potential 
^ at the origin. Thus the gauge potential has to be regularized. 

A possible regularization is provided by considering the vortex with the 
flux concentrated in a thin annulus of a radius R. Let us first consider the 
Hamiltonian in the range of parameters < MB < 1/4. The zero-energy 
state of the Hamiltonian outside the annulus is then a linear combination 



of the modes ^± given by Eq. (17). Inside the annulus the vector potential 



A = 0, and the zero-energy modes are 

*<r =Ci U M +C 2 i U } , (22) 

^Io(Kr) / I z/o(A-r) 



with A± given by Eq. (16), and C\ t 2 being complex constants. By matching 
these solutions at r = R, and taking R — > 0, we obtain, up to a normalization 
constant, the zero-energy state of the form 

p -X + r _ p -\-r I e -i<p 
nr) = ~ ^ . |. (23) 
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Notice that this zero-energy state is regular at the origin which is a conse- 
quence of the regularity at the origin of the solutions (J22]) of the vortex-free 
problem. Similarly, one may show that when MB > 1/4 the zero-energy 



mode is given by the spinor \l/ 2 i n Eq. (18) also regular at the origin and 
behaving ~ r 1 / 2 when r — > 0. 



6. Self-adjoint extension of the Hamiltonian 

Although the above regularization results in concrete solutions to the 
problem, we should consider the self-adjoint extension of the corresponding 
Hamiltonian Q in a more general manner by specifying the proper Hilbert 
space. That way, a family of Hermitian Hamiltonians is obtained, depending 
on free physical parameters that determine the scattering at the vortex core 
and the detailed profile of the single zero- mode (per spin). In the last section, 
we will comment on the regularization provided by the tight-binding version 
of the model, Eq. g. 

The application of the standard theory of self-adjoint extensions (SAE) 
[301 EZl ESI EH] prescribes that we need to ensure that the massive Dirac 
Hamiltonian (a differential operator) becomes Hermitian (self-adjoint) only 
after choosing the proper Hilbert space (i.e. domain of functions) on which 
it is allowed to act. Instead of analyzing the imaginary spectrum (which 
needs to be removed), we implement von Neumann's construction by looking 
directly at the conditions under which the Hamiltonian is Hermitian when 
acting on arbitrary functions that are square integrable (but might diverge at 
the origin, due to diverging potential there). This will effectively determine 
the coefficients of the linear combination C f i\t r ++C f 2'0- in angular momentum 



14 



channel I = and thereby fix the zero-mode. 

Using symmetry, we start from the radial part of the operator, H l (r), 
which acts in the subspace of angular momentum I spanned by functions of 
the form 

(e~ i<p ui \ 
I , (24) 

completely determined by ("J). Recall that the zero-energy states of the 
M — B model Q in presence of a 7r-flux vortex come in the form of Kramers 
pairs 





where ijj(r } tp) is exactly of the form in Eq. (24). 

We also implement the standard change of scalar product in r-space by 
rescaling ipi (r) = -^ipi(r), dipi{r) = -^dipiir) {d always denotes d/dr), where 
d = d — 7T-, after which H l takes the form 

1 ' ~ V -< d - l r) -M-B(02_i<i+i)) J ■ y M > 

We have reverted to the standard derivative (d) here. Since the 7r-flux enters 
through I — y l e ff, we omit it here. The non-derivative terms will not play any 
role in the following analysis since the standard 'centrifugal force' provided 
by I 7^ does not lead to singularities. However, the gauge potential will 
have the chance to provide us with the boundary condition exactly when 
I = 0. 

Now, for two arbitrary wavefunctions <p, tp which are determined by F = 
( f ) , U = ( " ) (we dropped index I), respectively, the condition of hermiticity 
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of H l (r) becomes (note the change in rdr): 



H-H* 



dr (j)(ryH l (r)ij(r) - / dr ip(r)*H l (r)(j)(r 



B \ f*d r u — d r f*u — g*d r v + d r g*v 
B 
0. 



i \ f*v + g*u 



F*a 3 d r U - d r F*a 3 U - iF*<j x XJ 



(0) 



(27) 



where the a Pauli matrices act on the two component functions, which vanish 
at infinity, and are evaluated at the origin (point r = 0) in the next-to-last 
line. 

There is a continuous family of restrictions on the behavior of square- 



integrable functions at the origin, such that (27) is satisfied, leading to the 



Hamiltonian which is Hermitian on such a chosen domain. The proper 
parametrization of the most general restriction on the allowed domains is 



achieved by using the linearity of (27). Namely, we define two linear opera- 



tors Ti, T2 which map arbitrary functions, i.e. the domain of H\ onto their 
value at the boundary, i.e. the space of complex two component vectors: 



r*:^(r)^(!gj 



(28) 



These operators are defined by (27) 



B 



F*a 3 d r U - d r F*a 3 U - iF*o x U (0) = (T 2 F, T 1 U) - (T 1 F, T 2 U). (29) 



Notice that this form can always be achieved due to the original form of the 
subtraction between H and H'. We can choose in particular, without loss of 
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generality, 

Tiff = Ba 3 d r U(0) - iytf(O), (30) 
T 2 U = U(0). (31) 

Any vector in the boundary space, i.e. ( £2 ) G C 2 , is an image by of some 
wavefunction, i.e. of some ( ~). Since U(0) = ff(ojl an d <9 r £7(0) = (^(oj) 
take on arbitrary values, this means that the boundary space indeed is = 
C 2 . 

The most general relation that has to be satisfied by a wavefunction such 



that (27) will hold is now parametrized by unitary mappings S in H b : 

V{Hu) ={tf>\{3- a )T^ + i(S + a Q )T 2 ^ = 0}, (32) 
with cr the 2x2 identity matrix, and T> denoting the domain of operator. One 



can directly understand from (29) that forcing arbitrary linear combinations 



of a general U(0) and d r ll(0) to zero will still preserve the condition (27), due 



to the linearity and the antisymmetric nature of the form of this expression. 



Eq. (32 ) is giving us a precise recipe and parametrization of the fact that this 
is the most general restriction that needs to be made on the wavefunctions 
U{r), i.e. on the ip(r). 

We now proceed to use the form of Tj to explore the allowed boundary 
conditions on the wavefunctions, in particular determining whether there is 
a self-adjoint extension Hu with the previously found zero energy states in 
its domain. For concreteness we focus on the case MB < 1/4. 

Since Tib is C 2 , our mappings 5 G (7(2), in contrast to the same problem 
in the case of an ordinary massive Dirac Hamiltonian necessitating a £7(1) 
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parametrization [31 J . The mappings S G U(2) can be parametrized by 



1 / e 1 U \ 



(33) 

with the quaternion basis = (o"o, ic?). We will use the notation [m , mi, m 2 , m 3 ] 
to represent the quaternion J^m^cr^, while to label the boundary states of 
the zero energy wavefunctions if) = Cx^ + + C^- we will use 

= *(O)=vWbU= I + ^ j, (34) 

\%{c x + c 2 ) y 

0' = ^~o(O) = a r (^o)| r=o = - ( ClXl + ° 2X2 ]. (35) 

\i(CiAi + C 2 A 2 )y 

Due to normalization, both and (f)' depend only on the vector Q), where 

X EE (36) 

It turns out that the case with 77 = is special, and so we examine it first in 
detail. Only this subclass contains the extension with regular zero-modes. 

6.1. Extensions described by SU{2), 77 = 



Eqs. (32), (33) lead to the following condition on two quaternions P, Q: 



B-P<P' = Q<f>, (37) 

mi 1 — itlq 7713 m 2 

Q = L-g - + H 1 + m o), — g h zmi ' "2" + Zm2 ' — 2~ + Zm ^' 

P = [—77713, —zm 2 , 77711, 7(7710 — 1)]. (38) 

Once the values of M, B (and therefore also Ai/ 2 = A±) are given, this equa- 
tion determines x as function of the particular SAE (if a solution for x 
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exists), and x then determines the specific linear combination of A1/2 decay- 
ing functions in the zero mode, influencing also its regularity properties at 
the origin. Since det(P) = 2(m — 1), we first consider: 

6.1.1. Extension with rj — 0, mo = 1 



According to definition (33) this immediately implies fa = 0, and 



0^ 



= 
Ci 



> \/r^o|r=o = 0, i.e. 

-c 2 . 



(39) 
(40) 



Such a wavef unction is regular at the origin and localized on the scale l/min(Aj), 
no matter the values M,B (or Ai^). We therefore see that such a wavefunc- 
tion is allowed when the physical SAE is given by 5* = [1, 0, 0, 0], 77 = 0. This 
will turn out to be the only extension allowing x = — 1, see Fig. [TJ essentially 
because it is the only extension for which the matrix on the left-hand side of 



Eq. (37) vanishes. 



6.1.2. Extension with 77 = 0, m < 1 



We must now consider Eq. (37) as a vector equation, treating Ci j2 as 
unknown variables. We get 



Y 



Ci 

Co 







Y 



Ba\\+p Ba\ 2 +p 
BbXx + q Bb\ 2 + q 

det(F) = B(X l -X 2 )(aq-bp), 



(41) 

(42) 
(43) 
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where a = P n + iP\ 2 , b = P 21 + iP 22 , P = Qii + iQi 2 , Q = Q21 + iQ 2 2- There 
are solutions only when 

aq = bp^m 2 = 0. (44) 

So in the case of m 7^ 1, the zero energy mode is allowed in principle only 
for Hamiltonians with m 2 = 0. The value of x in such a system is given by 

x = -4t4> 9W = BaX+p = BbX + q, (45) 

where 

g(X) = B(l-m -m 1 -im 3 )\+ Q - ij (m 1 -l-im 3 )+ Q + zj m (46) 

is a linear function of A. From this expression it seems that another way 
for x = — 1 (and the zero mode becoming regular at the origin for all model 
parameters M,B), is if mi = 1, where we assume Ai 7^ A2. However, this 
is not the case, because when ni\ = 1 only one component of the vector 



equation Eq. (41) can be satisfied, and not both component equations at the 
same time. 

Note that the extensions identified here by having m 2 = include the case 
from previous subsection, mo = 1. In fact, demanding that x — —1 directly 



from Eq. (41 ) implies that the columns of Y must be the same, which reduces 
to the demand a = b = =3- P = [0, 0, 0, 0] since Ai^ are non-degenerate. It 
is easy to see, using the definition Eq. ( [33] ) , that this is only possible in the 
above considered case m = 1. 
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6.2. Dependence on MB 

To analyze possible values for x for given model parameters M, B, we 
parametrize the SAE using hyperspherical coordinates 

= (cos ip, sin ip cos 9, sin ip sin 9 cos <p, sin ip sin 9 sin <£>) , (47) 

where the hypersphere SU(2) ~ S 3 is properly covered by ip G [0, 7r], G 



[0, 7r], and <p G [0, 27r). The condition = 0, Eq. (44), means (f = tt/2 



and also ip ^ (m 7^ 1), (ip, 6) ^ (f,0) (mj 7^ 1) in this subsection 



This parametrization can be plugged into the function g(X) from Eq. (46) 
together with the values of A1/2 = 1=l=v/ 2 1 ~ 4fc , where k = MB, so that one 
can deduce x from the model parameters M, B and the boundary condition 
determined by ip, 6. Analysis shows that x(ip, 6) weakly depends on the value 
of k = MB, and so Fig. [T] presents the typical behavior for fixed k. On the 
(ip, 9) plane, the function x(ip, 9) is strongly localized and a sharply peaked 
dipole. One can use the real and imaginary parts of x to identify structure 
in the scattering phase shift coming from the regularized vortex. 

6.3. From SU(2) to U(2) extensions, rj 7^ 



The introduction of 77 7^ in Eq. (33 ) does not change the analysis method 
of previous subsections. 



To start, the expressions for quaternions P, Q in Eq. (37) are slightly more 



complicated, so that demanding x = —1 through the vanishing of left-hand 



side of Eq. (37), as in Eq. (39), becomes 



det(P) = 0, i.e. (48) 

777,3 -rn ) 

exp(M7) = -. 

m 3 - i(l - mo) 
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However, even when the determinant vanishes, the matrix P can not become 
zero matrix, unless cos (77) = 1 =>■ r\ = 0, which actually just takes us back to 
a SAE considered in previous subsections. 



Next, considering Eqs. (37) by treating Ci, C2 as unknowns, as in Eq. (41 ) 



implies that zero-modes can exist only for a SAE that has 

det(Y) = => (49) 

exp(M?) = — — — 50 

1 + m\ + «m 2 

fri + n7i\ sm(ip) sin(0) cos(0) . . 

tan J = 1 + sinty) cos(0) ' (51) 



where in the last line we used the parametrization from Eq. (47). 



We emphasize that here, as for 77 = case, demanding a zero-mode with 



x = — 1 implies from Eq. (41) that the columns of Y have to be the same. 
Since Ai^ are not degenerate for physical values of MB, this gives a = b = 0, 
i.e. the already considered demand that the matrix P vanishes, which is 
impossible when 77 ^ 0. 



Following the derivation of Eq. (45), the x in general is again given by a 
ratio 

x = - fe^Ag) (52) 
9{m IJ ,,r],X 1 ) 

where we omit the lengthy explicit expression for g. Notice that even when 
this ratio approaches a singular limit by g — > or 00, the previous paragraph 
ensures that x will not reach the value —1. 

To summarize, the SAE analysis of the M — B model shows that a U (2) 
parametrization which completes this model allows the existence of zero mod- 
esfor a subclass of Hamiltonians described by three free angle parameters, 



e.g. T],0,ip in Eq. (47). These determine a matrix relating the value of the 
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K 



n 



Figure 1: Vortex singularity regularization. The value of x = C1/C2, where the zero- 
mode wavefunction ~ r~ 1 / 2 (xe~ XlT + e~ A2r ), is presented as a function of the SAE 
parametrized by two angles (ip,9), and r\ — (see Section [6]). The wavefunction becomes 
regular at the origin for x = — 1. Notice that Re(x) approaches —1 for -0 = 0, where 
also Im(i) = 0. The plots are evaluated for MB = 0.15, but are representative for all 
topologically non-trivial model parameters < MB < 1/4. 



wavefunction spinor to its derivative, both taken at the origin, according to 



Eq. (37). There is only a single SAE, Section 6.1.1 , determined by i] = ip = 



which allows for the existence of a zero energy state which is both localized 
and regular (vanishing) at the origin. 



7. Comparison to the tight-binding M — B model 

In this Section we present the results from a tight-binding M — B model, 
which in momentum space has the limit Eqs. Q, We numerically study 
this model on a 31x31 sized square lattice with periodic boundary conditions 
and a 7r-flux — anti-7r-fTux pair positioned at maximal distance. The details 
of this tight-binding model written in real space, i.e. on the lattice, are 
presented in Ref. [TT] . 

We first isolate the zero energy mode localized on a single 7r-flux, using 
the nearly degenerate inversion-symmetric and anti-symmetric zero modes on 
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-3-11 3 
Re(x) 



Figure 2: a) Comparison of the vortex-bound zero- mode localization lengths £± (black 



lines) predicted by the present M — B model (see text after Eq. (21)), and measured 
in the tight-binding version of the model on a 31x31 lattice. Inset shows in detail the 
excellent agreement in the MB < 1/4 regime. The lattice model has a phase transition 
at Mi att / Bi att — 4 which is absent in the continuum model, since the gap closes at finite 
momentum, b) The complex constant x = C\ jCi determines the form of the zero- mode 
wavefunction. When x = — 1 the mode is regular at the 7r-flux, and this situation is 
realized throughout the regime MB < 1/4 (dots in the legend mark the values of MB of 
the plotted points). 

the finite lattice. The lattice symmetries guarantee the (exp (—i(f)),i) T form 
of the wavefunction spinor, with the polar angle, just as in the continuum. 



By fitting the radial envelope of the spinor using Eq. (17) and the angle- 



averaged wavefunction, we extract the two localization lengths ^f att in units 
of the lattice constant. The results are shown in Fig. [2]^a). The agreement 



with the present continuum model prediction, see text after Eq. (18) and 



Eq. (16), is excellent in the MB < 1/4 regime (inset of Fig. ^ a)). For 



MB > 1/4 the oscillatory part of the wavefunction (see Eq. (21)) makes the 



numerical fitting less reliable, and the agreement is only qualitative. 
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One should note that the lattice tight-binding model has a natural length- 
scale, the lattice constant a, beside its dimensionless constants Mi att , Bi att (in 
units of inter-site hopping energy A), see Eq. 0. In the present continuum 
model, Eq. Q, the parameters M, B determine two length-scales, B and 
5 1 / 2 M~ 1,/2 , where the former acts as a length-scale in the zero-modes. The 
lattice constant is simply absorbed by Mi att , Bi att to give M, B of the con- 
tinuum model (see Eq. ([3])), and that is why the agreement in Fig. |2^a) is 
quantitatively precise. (In the lattice calculation we actually set Bi att = 1 
and vary only Mi aU , driving therefore both the parameter MB = Mi at tBi aU 
and the lattice topological transition parameter Mi att / Bi att .) 

To determine the SAE which is realized in the lattice model, we calculate 



the complex constant x = C1/C2, which was defined in Eq. (36), and rep- 
resents the ratio of contributions of two singular functions in the zero-mode 



wavefunction, Eqs. (17), (21). Fig. |2j(b) shows that throughout the regime 
MB < 1/4, x keeps near the value x = —1, which is the special case of zero- 
mode being regular at the 7r-flux position. As shown above, only a single SAE 
allows this, and this SAE is also realized in a thin-solenoid regularization of 
the continuum model. 



8. Quantum numbers of the zero-energy modes 

We will now show that these zero-energy modes carry non-trivial charge 
or spin quantum number depending on their occupation. For that purpose, 
we will write the continuum 4x4 Hamiltonian (|3| coupled to a U(l) vector 
potential as 

H = holiih + Ai) + (M — B(k + A) 2 ) 7o (53) 
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with the vector potential A given by Eq. ([5]). Note that the unitary matrices 
73 = 02 <S> T2 and 75 = 02 <8> T\ anticommute with the gamma-matrices 7 a , 
a = 0,1,2. Therefore, the Hamiltonian anticommutes with the matrices 
T 3 = 27073 and T 5 = 27075 which then generate chiral (spectral) symmetry 
relating states with positive and negative energies, i.e., if H\E) = E\E), 
then, for instance, T^\E) = | — E), and the matrix reduces in the zero- 
energy subspace of the Hamiltonian. These two properties then imply that 
in the ground-state the expectation value of a traceless Hermitian operator 
Q is given in terms of the zero-energy states of the Hamiltonian 

(Q} = 1[ £ - J2 *Q*) ■ ( 54 ) 

\occupied unoccupied / 

Therefore, in the case of a quantum spin Hall insulator threaded by a 7r-fiux 
vortex, depending on the occupation of a pair of zero-energy states, there are 
four possibilities for the ground state quantum numbers. Namely, when both 
states are occupied or empty, according to the above expression, the charge is 
+e or — e and the spin quantum number is zero. On the other hand, when one 
of the states is occupied, the spin quantum number is +1/2 or —1/2, while 
the charge is zero. In that way, the spin-charge separation, characteristic for 
one-dimensional systems [22], appears also in a two-dimensional system [38J, 
and is tied to a topologically non-trivial nature of the quantum spin Hall 
state P21 [13]. Similarly, zero-energy modes bound to the vortex core in an 
antiferromagnetic state on a honeycomb lattice lead to the phenomenon of 
spin fragmentation [39J; see Ref. [10] for a general discussion of this class of 
problems in the context of Dirac systems. 

Finally, let us note that the topological stability of 7r-fiux zero-modes, 
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i.e. stability under smooth deformations of the vector potential, has been 
strictly proved for both Dirac and Schrodinger Hamiltonians jT6l [TTj , but for 
the M — B model, which is a sum of both, we are not aware of an analogous 
proof based on an index theorem. 

9. Conclusions 

In conclusion, we have shown that a pair of zero-modes bound to a mag- 
netic 7r-flux are a generic feature of the M — B model in the topologically 
non-trivial phase. The continuum Hamiltonian of that model has both Dirac 
and Schrodinger kinetic terms and a mass, and has not been studied in de- 
tail previously. We have analytically found the zero-modes in the presence of 
the 7r-flux vortex in the entire range of parameters describing a topologically 
non-trivial phase with the bandgap opening at the zero momentum. These 
modes are exponentially localized around the vortex core, and a particular 
regularization of the vector potential corresponding to the vortex yields the 
modes regular at the origin, but in general, as we have shown, the form 
of the solution depends on the short-distance regularization of the vortex. 
Vortex zero modes obtained within a lattice tight-binding model match the 
ones found for a particular self-adjoint extension of the continuum Hamilto- 
nian. Finally, we discussed a realization of the two-dimensional spin-charge 
separation through the vortex zero-modes. 
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